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The formulation of the Einstein field equations admitting two Killing vec- 
tors in terms of harmonic mappings of Riemannian manifolds is a subject in 
which Charlie Misner has played a pioneering role. We shall consider the hy- 
perbolic case of the Einstein-Maxwell equations admitting two hypersurface 
orthogonal Killing vectors which physically describes the interaction of two 
electrovac plane waves. Following Penrose's discussion of the Cauchy prob- 
lem we shall present the initial data appropriate to this collision problem. We 
shall also present three different ways in which the Einstein-Maxwell equa- 
tions for colliding plane wave spacetimes can be recognized as a harmonic 
map. The goal is to cast the Einstein-Maxwell equations into a form adopted 
to the initial data for colliding impulsive gravitational and electromagnetic 
shock waves in such a way that a simple harmonic map will directly yield 
the metric and the Maxwell potential 1-form of physical interest. 

*for Charles W. Misner on his 60 th birthday 
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1 Introduction 



Charlie Misner was the first to recognize that the subject of harmonic map- 
pings of Riemannian manifolds finds an important application in general 
relativity. In a pioneering paper with Richard Matzner |l] he found that 
stationary, axially symmetric Einstein field equations can be formulated as 
a harmonic map. Eells and Sampson's theory of harmonic mappings of Rie- 
mannian manifolds [[| provides a geometrical framework for thinking of a set 
of pde's, in the same spirit as "mini-superspace" that Charlie was to introduce 
H for ode Einstein equations a little later. The subsequent developement of 
the subject of space-times admitting two Killing vectors, that eventually led 
to its recognition as a completely integrable system 0-0 has employed 
another formulation of the stationary, axi-symmetric field equations due to 
Ernst H which is equivalent to that of Matzner and Misner. 

Charlie's later work on harmonic maps H encompasses a scope much 
broader than this specific problem and its power and elegance is bound to 
make a major impact on theoretical physics. 

I was privileged to be in contact with Charlie's ideas at that time and 



worked on the two Killing vector problem [|HJ], ||TT| . It was the hyperbolic ver- 
sion of gravitational fields admitting two Killing vectors that attracted my at- 
tention. This is the problem of colliding impulsive plane gravitational waves 
for which Khan and Penrose had presented a famous solution fl2| . My work 



finally led to the exact solution for colliding impulsive plane gravitational 
waves with non-collinear polarizations [13|] which is physically the most gen- 
eral solution of this type. It turned out that the Matzner-Misner formalism 
was the one more readily amenable to the hyperbolic problem, even though 
it was originally intended for the elliptic case, whereas the Ernst formulation 
fitted the elliptic problem, ie the exterior field of rotating stars, best. The 
relationship between these two formalisms is given by a Neugebauer-Kramer 



involution WA 



A few years after the solution |rj[ appeared, there was a remarkable 
avalanche of papers on colliding plane gravitational waves. There were im- 
portant papers examining the singularity structure of spacetimes resulting 
from the collision of gravitational waves |18| - However, there was also 
a mass of new exact solutions which are all essentially devoid of any physical 
interest because their authors had chosen not to solve the Cauchy problem 
with the initial data appropriate to generic plane waves, but rather they 
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started with a "solution" and derived (!) the initial data. This type of de- 
rived initial data for the collision problem describes some very peculiar plane 
waves indeed. An inordinately large number of such references can be found 
in fl5j. 

Nevertheless, physically interesting colliding wave problems are still open 
and waiting for an exact solution! Remarkably enough, the interaction of 
plane impulsive gravitational and electromagnetic shock waves is in this cat- 
egory. We have the Khan-Penrose and Bell-Szekeres [ljj solutions describing 
the interaction of either two impulsive gravitational, or two electromagnetic 
shock waves alone and also the solution of Griffiths ITTIl for the interaction 



of an impulsive gravitational wave with an electromagnetic shock wave. But 
the generic case where we must consider the collision of both type of waves is 
missing even in the case of collinear polarization. The important open prob- 
lem here is the construction of an exact solution of the Einstein-Maxwell 
equations that reduces to all, the Khan-Penrose, Bell-Szekeres and Griffiths 
solutions. There are various unsatisfactory treatments of this problem in the 



literature pTJ, 122]. I shall give its harmonic map formulation. 



2 Initial Data 



The problem of colliding plane gravitational waves was proposed and in 
essence solved by Penrose |23| in 1965 even though most people writing on 
this subject do not seem to be familiar with it. We shall use Penrose's for- 



mulation of the Cauchy problem [24], [25] to discuss the interaction of two 



plane waves, every one of which will consist of a superposition of an impul- 
sive gravitational wave and an electromagnetic shock wave. The interaction 
will be determined by an integration of the Einstein-Maxwell equations with 
initial data defined on a pair of intersecting null characteristics. The initial 
values of the fields will be those appropriate to a plane wave which is given 



by the Brinkman metric p6 ], 



ds 2 = 2du'dv' - dx' 2 - dy' 2 + 2 H(v' ,x' ,y') dv 



I 2 



(1) 



and the superposition of an impulsive gravitational wave and an electro- 
magnetic shock wave, with amplitudes proportional to a, b respectively, is 
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obtained for 



H = % (V 2 - x' 2 ) «y(«0 - y (>' 2 + y' 2 ) (2) 



where 5 is the Dirac delta-function and 6 is the Heaviside unit step-function. 

The Brinkman coordinate system employed in eq. ([I]) is useful because the 
superposition of waves travelling in the same direction is obtained simply 
by addition. However, the Brinkman coordinates are not suitable for the 
collision problem because of the explicit dependence of the metric on x',y'. 
For this purpose we must transform to the Rosen form where the metric 
coefficients will depend on v alone. This is accomplished by the Khan-Penrose 
transformation 

v 1 = v , 

= u + \x 2 F F v + \y 2 G G v , 



u: 



r rF (3) 

y' = yG, 



which results in 



provided 



ds 2 = 2dudv - F 2 dx 2 - G 2 dy 2 , (4) 

F vv = [-a8(v)-b0(v)\ F , 

G vv = [-a8(v)+b9(v) \ G. 1 ' 

These are linear, distribution-valued ordinary differential equations which 
can be solved using the Laplace transform 

rod 

F(s) = / e sv F(v)dv (6) 

JO 

and from eq.([5]) we find 

Hs) = -^ V2 [(s~a)F(0) + F v (0)] (7) 

where of F(0) , F v (0) are initial values. They are obtained from the continuity 
of the metric and its first derivatives across v — which requires F(0) = 1, 
F v (0) = 0. In this case inverting the Laplace transform we get 

F = cos{bv9{v)) - | sin(bv9(v)) (8) 
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and the result for G is obtained by letting a — > —a in eq.(||) as indicated by 
eqs.(|). 

In Rosen coordinates the general form of the metric that admits two 
hypersurface orthogonal Killing vectors is given by 

ds 2 = 2e~ M dudv - e~ u (e v dx 2 + e~ v dy 2 ) (9) 

where U, V, M depend on only u, v and comparison with eqs.M) and (0) shows 
that for the initial value problem the data is given by 

2 

e~ u = cos 2 (bv9(v)) - sin 2 (bv9(v)) 

_ v b + a tan (buO(v)) (10) 
b — a tan (bv9(v)) 
e~ M = 1. 

The limiting values of this result are familiar. If we have just an impulsive 
gravitational wave, we must pass to the limit b — > which yields 

e~ u = l-av 2 9(v) 

e~ V = | + avG J} V \ (11) 
1 — av9{v) K ' 

e~ M = 1 

as in the case of Khan and Penrose. Furthermore, in the limit a^Owe have 
only an electromagnetic shock wave 

e~ u = 1 — sin 2 (bv9(v)) 

e~ v = 1 (12) 
e~ M = 1 

which is the result for the Bell-Szekeres case. 

Spacetimes describing colliding plane waves are divided into four regions: 

Region I : u < 0,v < empty space before the collision 
Region II: u > 0, v < a plane wave 

Region III: u < 0, v > another wave travelling in the opposite direction 
Region IV: u > 0, v > the interaction region 

The initial values given above are on v = 0, the boundary between Regions 
III, IV and similar results hold on u = 0, the boundary between Regions II, 
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IV, determining the u- dependence. In the latter case the amplitudes of these 
waves will be given by different constants, say a — > p and b — > q, cf eg. fl28|) 
in sequel. The case of Griffiths' solution is a mixture where we have eqs.flTT]) 
between Regions III, IV and eqs.(|l|) with u replacing v on the boundary 
between Regions II, IV. 



3 Einstein-Maxwell Equations 

The Einstein-Maxwell field equations governing the interaction of two plane 
waves is well-known |[2g|| . Starting with the metric @ and the Maxwell 
potential 1-form A 

A = Adx (13) 

where A depends only on u, v, we find a set of Einstein field equations which 
can be grouped into two categories. First we have the initial value equations 



2U VV -U V 2 + 2M V U V -V V 2 = 2Ke u ~ v A 2 
2U UU -U U 2 + 2M U U U -V U 2 = 2Ke u ~ v A 2 



(14) 



and their integrability conditions 



2A UV V V A U V U A V = 
2V UV - U V V U - U U V V + 2Ke u ~ v A u A v = 
2M UV + U V U U ~ V V V U + 2Ke u ~ v A u A v = 



(15) 



where k is Newton's constant in geometrical units. 

In eqs.(^) we have the wave equation for e~ u and its solution is im- 
mediate from the initial values. The following two equations are the main 
equations and the last equations is irrelevant as M can be obtained from 
quadratures once the main equations are solved. 

The problem consists of finding a solution to eqs. (|U|) satisfying the initial 
data (|T0|). Finally, we shall remark that eqs.flTTD, or (|12|) can be regarded 
as the solution of an initial value problem themselves, namely one between 
Region I and either a gravitational impulsive wave, or an electromagnetic 
shock wave across the null plane v = in Region III. In this mini-problem 
the first one of eqs. (|TJ|) serves as the field equation. 
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4 Harmonic Maps 



We refer to and [29] for a review and survey of the principal results on 
harmonic mappings of Riemannian manifolds. Here we shall briefly recall 
the most basic definitions in order to fix the notation. We shall consider two 
Riemannian manifolds endowed with metrics 

ds — Qziz d x d oc • % — X > > . . > t\i i , _ > 

ds' 2 = g' af3 d y a d y 13 , a = l,...,n' ^ ' 



and a map 



f:M^M' (17) 



between them. This map is called harmonic if it extremizes the energy fuc- 
tional of Eells and Sampson, 5X = 0, 



I ^ = J<^w^ urx (18) 

where the Lagrangian consists of the trace with respect to the metric g of 
the induced metric / * g on A4. When the target space M! is 1-dimensional, 
harmonic maps satisfy Laplace's equation on the background of M. and on 
the other hand if Ai is 1-dimensional, then harmonic maps coincide with the 
geodesies on M.' . The nonlinear sigma model corresponds to the harmonic 
map f :R 2 ^ S 2 . 

There are at least three different ways in which the Einstein-Maxwell 
equations (fl5|) can be formulated as a harmonic map. 

1. We can take a 4-dimensional target space M.' with the metric 

ds' 2 = e~ u hdMdU + dU 2 - dV 2 ) - 2^e~ v dA 2 (19) 

which has sections of constant curvature and the flat metric 

ds 2 = 2dudv (20) 

on M.. This is the electrovac analogue of the formulation given in [|D|]. It 



is not the most economical approach because M, which can be obtained by 
quadratures, appears explicitly in the metric (|19D. 

2. We can get rid of M and consider a 2-dimensional target space with 
constant curvature at the expense of regarding U as a given function on 
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A4 which does not enter into the variational problem as one of the local 
components of the harmonic map. Thus assuming that e~ u satisfies the 
wave equation, as in eqs.(|T5|), we can take the metric on the target space as 

8 = e -udSd£ 

\Re£f V ; 

where 

£ = e (V-U)/2 + l ^ A (22) 

is an Ernst potential type of complex coordinate. The metric (|21| ) is that of 
a space of constant negative curvature and £ is the complex coordinate for 
the Poincare upper half plane. There exists another representation, namely 
Klein's unit disk for the space of constant negative curvature. This is ob- 
tained by the transformation 

£ = J^L (23) 

and 

~ (i -CO 2 { } 

is the resulting form of the metric. Frequently this is the most convenient 
representation. The original Matzner-Misner ITT] as well as the Neugebauer- 



Kramer [|TJ]] formulations are of this type. 

The metric on Ai is the same as the one in eq. (pOf ). 

This approach is by far the most common procedure followed in the lit- 
erature. 

3. It is possible to reformulate the reduced problem by avoiding the ad 
hoc introduction of U on A4' at the expense of adding a new Killing direction 
to M. where the magnitude of the Killing vector is e~ u . In this case we have 

ds' 2 = , di %^ (25) 



where the definition of £ is the same as before and 

ds 2 = 2dudv - e~ 2U dz 2 (26) 
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where we consider the mapping to be independent of z and once again U 
is specified a priori. It appears that this possibility has not been discussed 
before in the literature. 

The first formulation where it is not necessary to introduce information 
from outside the variational principle is attractive. However, the close re- 
semblance of the latter formulations to non-linear a-models is an advantage. 



5 Solutions 

The metric (|20|) on M is not in a form most suitable for the construction of 
harmonic maps. We need to rewite it using new coordinates in such a way 
that some information about the initial data ( |10D is already incorporated 
into the system with the result that we can look for simple harmonic maps 
automatically satisfying the initial data. For this purpose we start with the 



wave equation for e and following Szekeres |28[ write its solution as 

e~ u = f(u)+g(v) (27) 
where from eq.(^) we know that 

7 / 2< (28) 

We can now consider a trivial coordinate transformation 

u -> f(u) v -> g(v) 

which amouts to a replacement of u, v by /, g in the Einstein-Maxwell equa- 
tions (r5|). We shall further introduce the following definitions which are 



formally the same as those given by Khan and Penrose 

(29) 



P = \l\~f Q = A-9 



\ + f w = ^\ + g 

in terms of which it will be convenient to introduce new coordinates on M. 
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In the vacuum case with non-collinear polarization we had found that 

a = pw — qr t — pw + qr (30) 

were useful new coordinates because of two reasons. First of all, a and r can 
be recognized as prolate spheroidal coordinates [[HJ and we have the Kerr- 
Tomimatsu-Sato solutions of the main equations. Furthermore, the simplest 
and the most familiar solution of this type 



£ = cos(a)r + isin(a)a 



(31) 



satisfies the initial data. 

This situation changes for the electrovac case. It turns out that a useful 
definition of new coordinates is 



a 



P 



p + q 



r + w r + w 

in terms of which the metric on M is given by 



ds : 



ti- 



de 2 



(1 + r 2 ) 2 

and from eqs.(|27|) and (|28|) we have 



1 + 



2\2 



(32) 



(33) 



. — p 2 — q 2 . 



(34) 



Then the main equations become 

r 2\2 



{l -a 2 ) m 
_ 2£ 



I -a 1 



Ai + o 2 ) 2 



+ 



1 + T 



2\2 



d 

:i-r 2 ) ~dr 



+ (1 + T 



2\2 



'1 T 2 ) 9 ^ 



aey 



(35) 



which is similar to the prolate spheroidal case but differs from it in some 
important respects. Its advantage lies in the fact that 



±1 



(36) 
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is a solution of eq. ([35]) that leads to the Bell-Szekeres solution. Furthermore 
it can be readily verified that £ = er is also a solution as in eg . (|3l|) which 
is again the Bell-Szekeres solution. In terms of these coordinates the Bell- 
Szekeres solution is given by 

(1 + r 2 ) 2 (1 + a 2 ) 2 \l+r*J \l + a 2 J y K J 

which may help to clarify the meaning of a and r. 

So the remaining problem is to a find a one-parameter complex solution 
of eq. (|35|) that reduces to eq.([3~6"l). Such a solution will be of physical interest 
as it will automatically satisfy the proper initial data. 

6 Conclusion 

The only proper conclusion of a paper such as this, namely the exact so- 
lution of eqs.(|35"D satisfying the initial data in eqs. fllOD is missing. In this 
paper I have given a list of the essential properties that we must require 
from a physically acceptable solution describing the interaction of plane im- 
pulsive gravitational and electromagnetic shock waves and presented some 
preparatory material towards such a solution. In my case this solution has 
been missing since 1978 and that is why I felt it inappropriate to publish 
the work reported here earlier. However, I now feel that the abundance of 
so many irrelevant "solutions" of this problem in the literature has made the 
presentation of the real problem imperative. 
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